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Abstract. A symmetry group method is used to obtain exact solutions for a semilinear 
radial heat equation in n > 1 dimensions with a general power nonlinearity. The method 
involves an ansatz technique to solve an equivalent first-order PDE system of similarity 
variables given by group foliations of this heat equation, using its admitted group of scaling 
symmetries. This technique yields explicit similarity solutions as well as other explicit 
solutions of a more general (non-similarity) form having interesting analytical behavior 
connected with blow up and dispersion. In contrast, standard similarity reduction of this 
heat equation gives a semilinear ODE that cannot be explicitly solved by familiar integration 
techniques such as point symmetry reduction or integrating factors. 



1. Introduction 

In the study of nonlinear partial differential equations (PDEs), similarity solutions are 
important for the understanding of asymptotic behaviour and attractors, critical dynamics, 
and blow-up behaviour. Such solutions are characterized by a scaling homogeneous form 
arising from invariance of a PDE under a point symmetry group of scaling transformations 
that act on the independent and dependent variables in the PDE [UI2]- 

For scaling invariant PDEs that have only two independent variables, similarity solutions 
satisfy an ordinary differential equation (ODE) formulated in terms of the invariants of the 
scaling transformations. However, this ODE can often be very difficult to solve explicitly, and 
as a consequence, special ansatzes or ad hoc techniques may be necessary in order to obtain 
any solutions in an explicit form. The same difficulties occur more generally in trying to 
find explicit group-invariant solutions to nonlinear PDEs with other types of point symmetry 
groups. 

An interesting example is the semilinear radial heat equation 

Ut = u rr + (n — l)r~ u r + k\u\ q u, k = ±1 (1) 

for u(t,r), which has the scaling symmetry group 

t A 2 t, r ->■ Ar, u -> \- 2/g u, (2) 

where r denotes the radial coordinate in n > 1 dimensions. This equation describes radial 
heat flow with a nonlinear heat source/sink term depending on a power q ^ 0. The coefficient 
k of this term determines the stability of solutions to the initial- value problem. In particular, 
for k = —1 all smooth solutions u(t,r) asymptotically approach a similarity form u = 
£ -1 ' 9 [/(r/\/t) exhibiting global dispersive behaviour u — > as t — > oo for all r > 0, while for 
k = 1 some solutions u(t, r) exhibit a blow-up behaviour u — > oo given by a similarity form 
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u — (T — t) l l q U (r/y/T — t) as t — ?■ T < oo [3] E]. In both cases, U satisfies a nonlinear 
ODE 



U" + ((n - 1)£ 



i 



\k£)U' - \qkU + kU\U\ q = 



(3) 



with 




r/y/i, k = 



r/y/T=t, k = 1 



-1 



(4) 



which arises from the scaling invariance (J2J). This ODE cannot be explicitly solved by 
standard integration techniques [2] such as symmetry reduction or integrating factors when 
and q ^ 0. (More specifically, ODE 03]) has no point symmetries X = rj(£, U)d/dU + 
({£,U)d/d£ and no quadratic first integrals = A{£,U) + B(£,U)U' + C7(£, [7)£/' 2 , as 
established by solving the standard determining equations [H [2] for X and \V.) As a result, 
few exact solutions U(£) other than the explicit constant solution U = (q/2) 1 / q are apparently 
known to-date. 

In this paper we will obtain explicit exact solutions for the heat equation (JTJ) by applying 
an alternative similarity method developed in previous work |5J on finding exact solutions 
to a semilinear radial wave equation with a power nonlinearity. The method uses the group 
foliation equations [6] associated with a given point symmetry of a nonlinear PDE. These 
equations consist of an equivalent first-order PDE system whose independent and depen- 
dent variables are respectively given by the invariants and differential invariants of the point 
symmetry transformation. In the case of a PDE with power nonlinearities, the form of the 
resulting group-foliation system allows explicit solutions to be found by a systematic sep- 
aration technique in terms of the group-invariant variables. Each solution of the system 
geometrically corresponds to an explicit one-parameter family of exact solutions of the orig- 
inal nonlinear PDE, such that the family is closed under the given symmetry group acting 
in the solution space of the PDE. 

In Sec. |2l we set up the group foliation system given by the scaling symmetry ((2]) for the 
heat equation (JTJ and explain the separation technique that we use to find explicit solutions 
of this system. The resulting exact solutions of the heat equation are summarized in Sec. [3J 
These solutions include explicit similarity solutions as well as other solutions whose form is 
not scaling homogeneous, and we discuss their analytical features of interest pertaining to 
blow-up and dispersion. Finally, we make some concluding remarks in Sec. HI 

Related work using a similar method applied to nonlinear diffusion equations appears in 
Ref . [7] . Group foliation equations were first used successfully in Refs. [8] [9] [10] for obtaining 
exact solutions to nonlinear PDEs by a different method that is applicable when the group 
of point symmetries of a given PDE is infinite-dimensional, compared to the example of a 
finite-dimensional symmetry group considered both in Ref. [5] and in the present work. 

2. Symmetries and group foliation 

For the purpose of symmetry analysis and finding exact solutions, it is easier to work with 
a slightly modified form of the heat equation ([T]): 




k = const., q ^ 0. 



(5) 
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In n > 1 dimensions, this heat equation (J5J) admits only two point symmetries: 

time translation X tra ns. = d/dt for all q, (6) 
scaling X sca i. = 2td/dt + rd/dr — (2/q)ud/du for all q ^ 0, (7) 

where X is the infinitesimal generator of a one-parameter group of point transformations 
acting on (t,r,u). There are no special powers or dimensions for which any extra point 
symmetries exist for equation (jSJ), as found by a direct analysis of the symmetry determining 
equations. 

To proceed with setting up the group foliation equations using the scaling point symmetry, 
we first write down the invariants and differential invariants determined by the generator 
(JTJ). The simplest invariants in terms of t, r, u are given by 

x = t/r 2 , v = u/r p , (8) 

satisfying X scal x = X scaX v = with 

V = -2/q. (9) 

A convenient choice of differential invariants satisfying x£^ G = X^ ff = for G(t,r,u t ) 
and H(t,r,u r ) consists of 

G = r 2 - p u t , H = r^V, (10) 

where X^ = X scal . — (2 + 2/q)u t d/du t — (1 + 2/q)u r d/du T is the first-order prolongation 
of the generator (JTJ). Here x and v are mutually independent, while G and H are related by 
equality of mixed r, t derivatives on u t and u r , which gives 

D r (r p - 2 G) = D t {r v - l H) (11) 

where D r ,D t denote total derivatives with respect to r, t. Furthermore, v,G,H are related 
through the heat equation ()5]) by 

r p - 2 G - Dr^H) = r p - 2 {(n - 1)H + kv q+1 ). (12) 

Now we put G = G(x,v), H = H(x,v) into equations (ITT]) and (fT2"|) and use equation (jSJ) 
combined with the chain rule to arrive at a first-order PDE system 

(p - 2)G - pvG v - 2xG x — H x + HG V - GH V = 0, (13) 
G — (p + n — 2)H + pvH v + 2xH x - HH V = kv q+ \ (14) 

with independent variables x, v, and dependent variables G, H. These PDEs are called the 
scaling- group resolving system for the heat equation (J5J). 

The respective solution spaces of equation (J5J) and system ( 1T5|) - (IT4^ are related by a group- 
invariant mapping that is defined through the invariants © and differential invariants (fTUj) . 
In particular, the map (G, H) — > u is given by integration of a consistent pair of parametric 
first-order ODEs 

ut = r p ~ 2 G(t/r 2 , u/r p ), u r = r p ~ x H{tjr 2 , u/r p ) (15) 

whose general solution will involve a single arbitrary constant. The inverse map u — > (G, H) 
can be derived in the same way as shown in Ref. jS] for the wave equation, which gives the 
following correspondence result. 
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Lemma 1. Solutions (G(x, v ), H(x, v )) of the scaling- group resolving system n~3\) -(T%\) are 
in one-to-one correspondence with one-parameter families of solutions u(t,r,c) of the heat 
equation |3J) satisfying the scaling-invariance property 

\- p u(\H, Ar, c) = u{t, r, c(A, c)) (16) 

where c(l, c) = c in terms of an arbitrary constant c. 

This correspondence leads to an explicit characterization of similarity solutions of the 
heat equation ([5]) in terms of a condition on solutions of the scaling-group resolving system 
( jT3]) -( TH]) . Consider any one-parameter family of solutions 

u(t,r) = r p v, v = V(x,c), (17) 

having a scaling-homogeneous form, where 

Ax V" - (1 + {2p + n - A)2x)V + p(p + n - 2)V + kV q+1 = (18) 

is the ODE given by reduction of PDE fl5]). From relation ffTUl) we have 

G(x,V(x,c)) = V'(x,c), H(x,V{x,c))=pV{x,c)-2xV'{x,c). (19) 

Next we eliminate c in terms of x and v by using the implicit function theorem on V(x, c)—v = 
to express c = C(x,v). Substitution of this expression into equation f[T9~j) yields 

H + 2xG = pv (20) 

where G = V'(x, C(x, v)), H = pv — 2xV'(x, C(x, v)) are some functions of x, v. The relation 
( 1201) is easily verified to satisfy PDE (EE3]). In addition, PDE P3J simplifies to 

- 4x 2 (G x + GG V ) + (1 + [2p + n - A)2x)G = p(p + n - 2)v + kv q+1 . (21) 

We then see that the characteristic ODEs for solving this first-order PDE are precisely 

dv/dx = G, -Ax 2 dG/dx + (l + (2p + n-A)2x)G = p(p + n - 2)v + kv q+1 , (22) 

which are satisfied due to equations ( |T8|) and ( TL9l . Hence, we have established the following 
result. 

Lemma 2. There is a one-to-one correspondence between one-parameter families of simi- 
larity solutions ( (77| ) of heat equation |3J) and solutions of the scaling-group resolving system 
/ TJ^) -p^]) that satisfy the similarity relation 

We now note that, under the mapping ( Tl5l) . static solutions u{r) of the heat equation 
correspond to solutions of the scaling-group resolving system with G — 0. Consequently, 
hereafter we will be interested only in solutions such that G^O, corresponding to dynamical 
solutions of the heat equation. 

To find explicit solutions of the PDE system ffT3l) - fTHl) for G(x, v), H(x, v), we will exploit 
its following general features. First, the power nonlinearity ku q+l in the heat equation 
appears only as an inhomogeneous term kv q+1 in the PDE fTH|) . Second, in both PDEs 
ffT3]) and fTH|) the linear terms that involve v derivatives have the scaling homogeneous form 
vG v and vH v with respect to v. Third, the nonlinear terms in the homogeneous PDE (fT51) 
have the skew-symmetric form HG V — GH V , while HH V is the only nonlinear term appearing 
in the non-homogeneous PDE (TH1) . These features suggest that this PDE system can be 
expected to have solutions given by the separable power form 

G = gi(x)v a + g 2 (x)v, H = /ii(x)v a + h 2 (x)v, a ^ 1. (23) 
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For such an ansatz, we readily see that the linear derivative terms G x , H x , vG v , vH v in 
each PDE (|T3|) and (I14p will contain the same powers v,v a that appear in both G and H, 
and moreover the nonlinear term HG V — GH V in the homogeneous PDE (TT3l will produce 
only the power v a due to the identities v a {v ) v — v(v a ) v = (a — l)v a and v(v) v — v(v) v = 
v a (v a ) v — v a {v a ) v = 0. Similarly we see that the nonlinear term HH V in the non-homogeneous 
PDE (JT3D will only yield the powers v, v a , v 2a ~ l . Since we have a / 1 and 5 ^ 0, the 
inhomogeneous term kv g+l must therefore balance one of the powers v 2a ~ x or v a . 

In the case when we balance q + 1 = a, the terms containing v a = v q+1 and v 2a ~ 1 in PDE 
f fT4|) immediately yield 

hi = 0, 9l = k. (24) 
Then the terms containing v a in PDE (fl3|) reduce to 

/i 2 = (25) 
which leads to a simplification of the remaining terms in both PDEs, yielding 

xg' 2 = g 2 = 0. (26) 



Thus, for this case, the ansatz ([23]) gives a one-term solution for G with H = 0. 
In the other case, balancing 2a — 1 = q + 1, we get 

a = l + §. (27) 
The terms containing t> , v a , ti 2a_1 = in the PDEs ( fl3l) and ( fl4l) then yield 

/ii + 2x^ + (1 - hi 2 ) 9l + |/i l52 = (28a) 

h' 2 + 2xg' 2 + 2g 2 = (28b) 

2xh[ -(n-l + ^-^h 2 )h 1 + g x = (28c) 
2xh' 2 - h 2 2 - (n-2)h 2 + g 2 = (28d) 
(l + |)/ il 2 + A; = 0. 
Through equations (I28e[) . f l28cj) and (128dj) . we obtain 



1 + |)/n 2 + fc = 0. (28e) 



gi = ± x l-^(n-l+ q -±±h 2 ) (30) 

^ 2 = -2x/i' 2 + (n - 2 + h 2 )h 2 (31) 

and thereby we find that equations (I28al) and (128bl) reduce to an overdetermined system of 
nonlinear ODEs 

Ax 2 h' 2 - 2xh 2 - (1 + 2(n - 2)x)/i 2 = c = const. (32a) 
4x/i 2 - g ^ 2 + ^ /^ 2 2 + 2/i 2 + n - 1 = (32b) 



where c is an arbitrary integration constant. The ODE system (|32|) can be solved by a sys- 
tematic integrability analysis, which we have carried out using computer algebra (discussed 
in more detail in Sec. 12. ip . The results of the analysis give three two-term solutions with 
G^O. 

Proposition 1. For n ^ 1, ansatz ( dSP yields altogether four solutions of the scaling- group 
resolving system $T3\I-[I$ with G / 0: 

G = kv q+1 , H = 0, (33) 

G = ±(4 - n) / -Mn-2) ^ (w _ 3)/(n _ 2) ^ = + _ , (34) 

V n — 3 A — n 

q = -^—^-\, 2,3,4 
2 — n 

G=^^±^, H=-G--, q = -A, n = 5/2 (35) 

3x + 1 2w ' 3 2' y ' 7 v ; 

3v (1 ± ^/^2kv) 1 3x - 1 

3x + 1 o 3x + 1 



Acme o/ these solutions satisfy the similarity relation l[20\) . 

Motivated by the success of the ansatz (1231 . we now consider a more general three-term 
ansatz 

G = gi(x)v a + g 2 (x)v b + g 3 (x)v, H = /ii(x)w a + h 2 (x)v b + h 3 (x)v (37) 

where 

a^b, a 1, b^l. (38) 

This ansatz leads to a more complicated analysis compared to the previous two-term ansatz. 
Specifically, the homogeneous PDE ( TT3]) now contains the power v a+b ~ l in addition to v, v a , v b , 
while the non-homogeneous PDE (fl4l) contains the further powers v 2a ^ 1 ,v 2b ^ 1 ,v q+1 . We 
determine the exponents in these powers by a systematic examination of all possible balances. 

Firstly, since q ^ in PDE ( I14p . v q+1 must balance one of v a ,v 2a ~ 1 ,v a+b ^ 1 . (Note, by 
the symmetry a <H- b in the ansatz (1371) . the other possibilities v b ) v 2b ^ 1 for balancing v q+1 
are redundant.) Secondly, t> 2a_1 can balance only t> 9+1 or v b due to conditions ( 1381) . and 
otherwise if t> 2a_1 is unbalanced then its coefficient ah\ 2 must vanish. Likewise f 26 " 1 can 
balance only v q+1 or v a , and otherwise its coefficient bh2 2 must vanish. In a similar way, 
either y a+b ~ l balances v q+1 or v , and otherwise if unbalanced then its coefficient 

(a + b)h\h 2 vanishes. Finally, in PDE ffl3|) . t> a+b_1 can balance only v, and otherwise the 
factor g\h 2 — g 2 hi must vanish in the coefficient of t> a+6_1 . 

Several cases arise from examining all of these different possibilities. After eliminating all 
trivial cases that lead to h 2 = g 2 = (whereby the ansatz (137]) just reduces to the previously 
considered two-term case (|2"3"|) ). we find the following non-trivial cases to consider: 

q = a = 2, 6 = 0; (39) 
g = -3/2, a = 0, b =-1/2; (40) 
q = -2/3, a = -b = -1/3. (41) 
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For case the PDEs flB]) and QH) yield 

/4 + 2^; + (i - /i 3 )^i + Ma = o 

2x/i; - (n - 1 + 3/i 3 )/ii + #i = 
h' 3 + 2xg' 3 + 2g 2 h 1 - 2 9l h 2 + 2g 3 = 
2xh' 3 - /i 3 2 - 2h 1 h 2 + (2 - n)h 3 + # 3 = 

/l' 2 + 2x^2 + (3 + ft. 3 )# 2 - #3^2 = 

2xh' 2 - h 2 h 3 + (3 - n)h 2 + #2 = 
2h 1 2 + k = 0. 
From equations fll2gD, fl42fj) . ( 142dp . fl42bj> . we have 



#2 
#3 



-4 

—2xh' 2 + (n — 3 + 7i 3 )/i2 

-2x/i 3 + (n - 2 + fi 3 )/i 3 ± V^2kh 2 



91 = ± 

and then equation (142 al) gives 

h 2 = < 



-k 



[n 



1 + 3/i 3 ) 



(4x/i 3 -2/i3 2 + 2/i3 + n-i; 



The remaining equations (142c[) and ( I42ej) become, respectively, 



4x 2 /i , 3 ' 



(12xA 3 + 2(n - 6)x + l)/i 3 + 4/i 3 3 - 6/i 3 2 + 2(3 - 2n)h 3 = 



(42a) 
(42b) 
(42c) 
(42d) 
(42e) 
(42f) 
(42g) 



(43) 
(44) 
(45) 

(46) 



(47) 



(48a) 



Ax^h'z" - x(Axh 3 + 2(n - 13)x + l)h' 3 ' + (Qxh 3 2 + (2(n - 12) x + l)/i 3 + 9(4 - n)x - 3/2)/i' 3 



- I2x 2 h' 2 - \(2h 3 2 - 2h 3 + 1 - n)(/i 3 2 - 2/i 3 + 5 - 2n) = 



(48b) 



which is an overdetermined system of two nonlinear ODEs for h 3 (x). We solve this system 
(148]) by an integrability analysis using computer algebra. This yields one solution with G = 0, 
plus two solutions with G ^ which are summarized in Proposition 2. 

In a similar way, each of the cases (HOI and (|4T|) leads to an overdetermined system of four 
nonlinear ODEs for h 2 (x), h 3 (x). For both cases the results of an integrability analysis yield 
only solutions with G — 0. 

Thus we have the following result. 



Proposition 2. For 1, ansatz (Sty yields two additional solutions of the scaling- group 
resolving system (TT3|)-([7^]] with G 7^ 0: 

2 



G = ±-V^2k v ± -== 
4 V 

15 ,—rf 1 



-G 
3 



G = ±—V^2k [vt 



H 



15 



G + w + 



2 
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Neither of these solutions satisfies the similarity relation §2U$ . 
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n = 5/2 
n = 5/2 



(49) 
(50) 



2.1. Computational remarks. 

The integrability analysis of the previous ODE systems is non-trivial due to the degree of 
nonlinearity of the ODEs and the algebraic complexity of the coefficients in addition to the 
appearance of parameters in each system. 

For the first ODE system ( I32p . the integrability analysis consists of the following main 
steps. We eliminate h 2 to get a single algebraic equation, which is quadratic in /i 2 . By 
differentiating this equation and using it to eliminate h 2 from either of the original ODEs 
in the system, we obtain a second algebraic equation, which is cubic in h^. The coefficients 
in each algebraic equation are expressions in terms of the independent variable x and the 
parameters q, n. We next use cross-multiplication repeatedly to eliminate the highest-degree 
monomial terms in both of the algebraic equations until one equation no longer contains \\2 
while the other equation is linear in h 2 . At each algebraic elimination step, we must note that 
a case distinction will arise if the coefficient of a highest- degree monomial vanishes for some 
values of q or n. For each case, once the final algebraic equations have been obtained, we 
solve the equation without h 2 by splitting it with respect to x, which will yield conditions on 
the parameters q, n, and we then solve the linear equation for h 2 subject to these conditions 
(if any). 

The integrability analysis for the second ODE system ([4*8]) is the same except that we 
must first use differentiation combined with cross- multiplication to eliminate h'£ and h'l, 
thus reducing the differential order of the system down to first-order, where the coefficients 
are expressions in terms of the independent variable x and the single parameter n. We 
may then proceed as before by using algebraic elimination to reduce this system to a linear 
equation that can be solved for /i 3 and an equation that does not contain h 3 and thereby 
determines n. 

A similar integrability analysis applies to the two other ODE systems, each of which 
requires solving four ODEs that contain two dependent variables I12, h% and their derivatives 
up to second order, in addition to the independent variable x and the single parameter n. 

Because of the complexity of the algebraic expressions and the number of case distinctions 
that arise in these analyses, it is very difficult for an automatic computer algebra program 
to fully classify and find all solutions. (For example, the Maple program RiffSimp running 
on a workstation for several days was unable to complete the full computation for any of the 
second, third, and fourth systems.) 

To overcome these difficulties, we have used the interactive package Crack [IT] which has 
a wide repertoire of techniques available, including eliminations, substitutions, integrations, 
length-shortening of equations, and factorizations, among others. Using Crack, the com- 
plete solution of the integrability analysis was obtained in about 50 interactive steps taking 
2 seconds in total for the system (1321) . and about 400 interactive steps taking 7 minutes 
in total for the system (jUJ), while less than 100 steps taking under 1 minute in total were 
needed for each of the other two systems [12]. 

3. Exact solutions 

To obtain explicit solutions u(t, r) of the heat equation ([5]) from solutions (G(x, v), H(x, v)) 
of its scaling-group resolving system (TlT^ - CHj) . we integrate the corresponding pair of para- 
metric first-order ODEs (1101) . The integration yields a one-parameter solution family u(t, r, c) 
which is closed under the action of the group of scaling transformations (T2]). 
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Theorem 1. The semilinear heat equation (TJJj has the following exact solutions arising from 
the explicit solutions of its scaling- group resolving system found in Propositions 1 and 2: 

u = (-kq(t + c))- 1/q , q^O (51) 



n-2 



(52) 

m= ^±v^ (1 + c(3t + r 2 )) (7 + r JJ ' <? = - 4 > " = 5/2 (53) 
5(3t + r 2 ) 

(r(15t + r 2 ) + Cv ^)^2fc 7 V ' 
3(t + c-r 2 ) 

u = ± i q = 2, n = 5/2 55 

r(3(t + c) + r 2 ) v /z 2fc ' K J 

5(3(t + c)+r 2 ) . . . 

u = ± V V ; I , 9 = 2, n = 5/2 56 

r(15(t + c)+r 2 ) v /z 2fc 7 V ; 

where c is an arbitrary constant. 

Modulo time-translations t — > t — c, solutions ( 15TT) . ( 1521 . ( 1551 . ( 156]) are similarity solutions 
since their form with parameter c = is preserved under scaling transformations (J2]) on 
r,t,u. In contrast, the parameter c in solutions (15"3"|) and (I54p cannot be removed by time- 
translations, and consequently the form of these solutions is not scaling-homogeneous since 
c gets scaled under the transformations ([2]) on r, t, u. Thus, with respect to the action 
of the full group of point symmetries generated by time-translations and scalings for the 
heat equation fl5]), the solutions ( 153]) and (J5JJ for c ^ yield two-dimensional orbits of 
non-similarity solutions given by 

u= (±^k(i + c (3(t + c)+r 2 )) ( 3(t ^ C) +rX\ ' , q = -4, n = 5/2 (57) 

5(3(t + c) + r 2 ) . . . 

u = ± V V ; ; ==, q = 2, n = 5 2 58 

(r(15(t + 5) + r 2 ) + c^) V^2k ' K J 



whereas the solutions ( l5Tj) . ( 1521 . (tool . (1561 represent one-dimensional time-translation orbits 
of similarity solutions 

u = {-kq)- l l q r- 2 / q {t/r 2 )- l l q 1 q^O (59) 
(n-2)(n-3)Y-" /2 , ^ n , ^,,--2 „ 2 



( ( " ^ 3) ) r"- 2 (±(l-2(n-4)t/r 2 ))"- 2 , 9 = ^- ± -1, n^2,3,4 

(60) 



M = ±(-2fc)- 1 /V- i y g = 2, n = 5/2 (61) 

3r/r z + 1 

M= ±(-2fe)-VV ^g±ll , , = 2, ,, = 5/2. (62) 



Note that solution ( 156|) is a special case of solution ( 158]) given by c = 0. We remark that 
solution ( 158]) previously has been obtained in work [13] on nonlinear diffusion equations 
through use of Bluman and Cole's nonclassical method. 

Among all these solutions, the ones ( 15 ip and (15 2 p that exist for integer values of n describe 
radial heat flow in W 1 . In section Sec. 13.11 we will discuss their analytical features related to 
blow-up and dispersion for n ^ 1. The remaining solutions ( 15"5"|) . ( ]5T)]) . ( ]5"Tj) . ( 15"5|) that exist 
only for a non- integer value n = 5/2 have a different interpretation describing heat flow in 
the plane M 2 with a point-source of radial heat flux at the origin, as we will show later in 
section Sec. 13.21 We will also discuss this interpretation for the solution fl52|) in the case of 
non-integer values of n. 

Before proceeding, we observe that the heat equation ([5]) for all values of n can be written 
in the form of a gradient flow 

Ut = -SE/Su (63) 

using the "energy" integral 

E = jf Qw 2 - kf(u)J r n ^dr (64) 

with 

/(«) = <? + 2 (65) 
I In q = —2. 

Here S/Su denotes the usual variational derivative (i.e. Euler operator) with respect to u. 
This integral (I64p obeys the equation 

J rp poo 

— = - lim(r n - V«t) - / u 2 t r n ~ l dr (66) 

J 

for any solution u(t,r) with sufficient asymptotic decay for large r. If the "energy flux" 
\mi{r n ~ l u r Ut) of a solution u(t,r) is non-negative then equation (I66p shows that E is a 

r— s>0 

decreasing function of t. As a result, for solutions u(t,r) that also are non-negative, both 
terms in E will be non-negative if k < and g 7^ —2, so then E decreases to zero as t — > 00. 
In this case u(t,r) will have dispersive behaviour such that u — > and w r — > for all r > 
as t 00. If instead k > or q = —2, then the two terms in E will have opposite signs, in 
which case E may decrease without bound, allowing u(t,r) to have blow-up behaviour such 
that u — > 00 or u r — > 00 as t — > T for some T < 00. 

3.1. Behavior of solutions for n = 2,3, etc. 

Similarity solution (15 ip is spatially homogeneous and has no restriction on the sign of k. 
It thus represents the general solution of the ODE u t = ku q+l with an arbitrary nonlinearity 
power q + 1 7^ 1. For q > the behaviour of u is determined by the sign of k. In the case 
k < 0, u = (\k\q)~ l l q / (t + c] x l q — > is dispersive as t — > 00, whereas in the case k > 0, 
it = (kq)-V 9 /(T - tf/ q ->■ 00 has a blow-up for i -> T with < T = — c < 00. 

In contrast, similarity solution (|52|) is restricted to the special nonlinearity power g + 1 = 
(n — 4)(n — 2) > 1/3 and requires k < and n > 5. In this solution, for all t > 0, u is 
singular at r = and is unbounded as r — >■ 00. Interestingly, u vanishes on the space-time 
parabola given by t + c = r 2 / (2n — 8) > 0. As a consequence we can modify the solution to 
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have better analytical behaviour by using this parabola (with parameter c = 0) as a cutoff 
such that 

(+ \ ) ' r - ftV ^ (R>7\ 

U (t,r) = \ / /_^n-2 n ( 67 ) 



(see figure [Q) with 



a 



p r 2-n ( t _ (r/ a ) 2 ) n , < r < «V« 
(n-2)(n-3)\ W2 



^2(n-4), /3 



|A;|(n-4) 2 



(68) 




Figure 1. Space-time graph for solution ([6 

At all space-time points away from r = 0, this modified similarity solution ( |6"7|) is con- 
tinuous and has continuous partial derivatives of order n — 3 > 2, so thus it satisfies the 
heat equation (jSJ) in the classical sense (i.e. it is at least C l in t and C 2 in r) on the spatial 



domain r > 0. At r = 0, u remains singular for all t > 0, but since r' 



n-l 



is non-singular as r — > 0, it is ra-dimensionally integrable on the whole spatial domain r > 
(see figure EJ). 

A physical interpretation of similarity solution ( 16?]) can be seen by considering the radial 
heat flux equation 

dH „ 

— = S + F (69) 



satisfied by the radial heat integral 



dt 



H 



ur 



n-l 



dr 



(70) 



n 



0.2 0.4 0.6 0.8 1 

r 

Figure 2. Graph of solution ([67]) with n = 6 (q + 1 = 1/2). 

which is a measure of the total amount of heat in u(t, r), where 

F = lim(-r n - V) (71) 

defines the outward radial heat flux at the origin, and 

POO 

S = k u q + x r n - x dr (72) 
Jo 

gives the net amount of heating or cooling produced by the nonlinear source/sink term in 
the heat equation (jSj). For the solution flBTl) these quantities are given by 

[" % (t - (rla) 2 ) n ' 2 rdr = — -Pt n ~ l = H > 0, (73) 
Jo n-1 

lim f-(2 - n)(3 (t - (r/a) 2 )™" 2 + 0{r)) =(n- 2)(3t n ~ 2 = F > 0, (74) 

[ OVi (3(^)/(n-2) ( t _ ( r / a fy-\*dr = -2(3t n ~ 2 = S<0. (75) 
Jo 

Hence u has a positive amount of heat (173|) that increases with t due to an increasing, positive 
radial outward heat flux at r = 0. Since n > 5, this flux (JTljl is greater than the net cooling 
(175]) caused by the nonlinear sink term. Therefore this similarity solution ( 16 7p physically 
describes the dispersion of heat produced by an outward radial heat flux at the origin in M. n , 
with some heat absorbed by a nonlinear heat sink proportional to -u*™ -4 )/*™" 2 ) at all points in 
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IR n . The dispersion has the behaviour of a radial temperature front at r = ayi that moves 
outward with speed dr/dt = a/(2\rt) for all t > 0. 

Interestingly, for long times, u increases without bound in t at any spatial point r > 
inside the temperature front. This non-dispersive temporal behaviour is related to u having 
an infinite "energy" ( 16%|) . i.e. E = +00, so that the flux equation (l66i) is not well defined, 
allowing u — > 00 as t — > 00 despite the "energy" integral being formally positive due to 
k < 0. 



3.2. Behavior of solutions for n / 2,3, etc. 

The heat equation §5$) can be written in a different form 

u t = u rr + (1 — v)r~ l u r + ku q+1 , v = const., k = const. (76) 

in terms of a parameter v = 2 — n which applies to non-integer values of n. This equation 
fT75|) describes radial heat flow in IR 2 with an extra source/sink term given by vu r jr [14J. To 
interpret this term physically, we consider the 2-dimensional radial heat integral 



oc 



H — urdr (77) 
Jo 

satisfying the radial flux equation 

^L = S + F + v\imu (78) 

dt r^O 

where 

F = lim(-rw r ) (79) 
defines the outward radial heat flux at the origin, and 

/■oo 

S = k u q+1 rdr (80) 



Jo 

gives the net amount of heating or cooling caused by the nonlinear source/sink term in the 
heat equation ( ITS"]) . The flux equation ( 1751) shows that, for the respective cases v > or 
v < 0, the term —vu r /r has the interpretation of a heating or cooling point-source at the 
origin in R 2 . Thus, solutions u(t,r) will physically describe radial heat flow arising from a 
point source in a thin layer, with heat also produced or absorbed at all points in the layer 
due to a nonlinear source/sink term ku q+1 . 

Note that the 2-dimensional heat equation (!76|) retains the form of a gradient flow (|63|) 
with n = 2 — v. 

Consider similarity solution ( 15 2 p with n = 2 — v: 

u(t,r) = (3r u (±(a(t + c) + r 2 ))~\ q = -2/v, v± 0,-1, -2 (81) 

where 

a = 2(i/ + 2), = + l)/kf/\ (82) 

which requires /c>0if0>z/>— 1 and k < if v < — 1 or v > 0. The behaviour of this 
solution depends essentially on the separate signs of a and v. 
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For a < 0, (IHTI) can be modified similarly to (|67|) by putting a cutoff on u at the space-time 
parabola where u = 0. Then 




•^h:,,,.,). ol^7R (83) 



gives a classical solution of the heat equation ( 1761) on the spatial domain r > (i.e. ii belongs 
to C 1 (IR + ) in t and C 2 (IR + ) in r), with v < —2 and A; < 0. However, at r = 0, u is singular 
such that, for all t > 0, H — +oo and F + u\mvu = 2{yH)t~ v ~ 1 lim r u+1 = +oo since 

r— >0 r—tO 

v + 1 < — f. The modified similarity solution ( )83|) therefore has the physical interpretation 
of heat dispersion produced by an infinite net outward radial heat flux at the origin in IR 2 , 
with a radial temperature front located at r = yjt/\a\ for t > 0. 

For a > 0, ( )8T|) is smooth and positive on the spatial domain r > 0. It thus gives a 
C°°(R x R + /{0}) solution of the heat equation (ITS"]) with the asymptotic behaviour 

/ x n u, 2x v ir u ((3a~ u t- u + 0(r)), r^0 
M (t,r) = /3r>t + r 2 r =<^ XrmH^/^ _w (84) 

Ir v (p + 0(l/r)), r — > oo 

(with c = 0) where z/ > —2. In the least interesting case when > v > —2, u is unbounded 
for large r, whereby it has an infinite amount of heat H = +oo and "energy" E = +oo. In 
contrast, when z/ > 0, u decays to for large r and vanishes at r = for t > 0. As a result, 
in this case the "energy" of u for t > is given by 

E = ^L_*V + 2u + 6){u + 2 )( 1 -)/ 2 r ( 1 +^^/2) /3t _ (1+3i/)/2 (g5) 

which is finite and positive due to A; < for v > 0. In particular, E decreases to as t — > oo 
in accordance with the flux equation f )66|) . Moreover, provided v > 2, u has a finite amount 
of heat 

H = ^2L {v + <X)^l\v - 2)' 1 r( / + Z/ wV W2 (86) 

2=w 2 3 v ; v ; r((i + z/)/2y v ; 

which decreases to as t — > oo. For t > the corresponding heat flux quantities in this case 
are given by 

S = - V -^H/t < 0, (87) 
F = 0, z/limu = 0. (88) 

r— >0 

Therefore, in the most physically interesting case v > 2, the similarity solution fl84|) 
describes an initial monopole-like circular heat distribution in M. 2 (see figure [3]) producing a 
smooth positive dispersive radial heat flow with some heat absorbed by a nonlinear heat sink 
proportional to v}~ 2 l v at all points in M 2 . In particular, this solution has no point-source or 
heat flux at the origin for all times t > 0. 

Similarity solutions flo3]) and flo^j) have q + 1 = 3, v — —1/2, and k < 0. In both 
solutions, for all t, w is singular at r = and has slow radial decay such that u = 0(l/r) 
for large r. Consequently, the amount of heat in u is H = ±oo. In particular, at t — 0, 
m is given by a 1/r heat distribution for r > 0, while for £ — >■ oo at any spatial point 
r > 0, u = r _1 (±(l/v / — 2k) + 0(l/t)) approaches a constant multiple of the same 1/r 
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Figure 3. Graph of solution (jHlj) with v — 3. 



heat distribution. This non-dispersive temporal behaviour occurs because u has "energy" 
E = +oo, which does not decrease with t. 

Non-similarity solution ( 154"1) also has q + 1 = 3, v — — 1/2, and k < 0. Compared with 
the previous similarity solutions fl55|) and f )56|) . it exhibits the same long-time behaviour 
it = r~ 1 / \/ — 2k) + 0(l/t)) as t — > oo at any spatial point r > 0. It also exhibits the same 
radial decay u = 0(l/r) as r — > oo for all t, whereby = ±oo and E = +oo. However, near 
r = 0, dMD is less singular than ([55]) and ([56]), such that u = r~ 1 / 2 (±(15/c v / -2A;)t + O(r)). 
As a consequence the heat flux quantities for ( 1541) have the properties F + v lim r _>o u = and 
< — S < oo. Moreover, at t = 0, (154"]) reduces to the heat distribution u = ±(5/ \J—2k){r + 
c/r 3 ' 2 ) -1 which is C 1 for r > if c > and which vanishes at both r = and r = oo 
(see figure H]). This non-similarity solution ( l54j) therefore has the physical interpretation of 
an initial circularly peaked heat distribution which vanishes both at the origin and spatial 
infinity in M 2 , producing a radial heat flow that is singular at the origin and has a 1/r decay 
for large radius, with a finite amount of heat absorbed by a nonlinear heat sink proportional 
to u 3 at all points in M 2 . The singularity in the heat distribution u at the origin corresponds 
to an infinite cooling point-source plus an infinite outward heat flux, whose net contribution 
to the heat flow for all times t > is zero. 

Finally, non-similarity solution ( 15 7p has q + 1 = —3, v = —1/2, and k > 0. In this solution, 
for all t, u is unbounded as r — > oo and has a cusp at the space-time parabolas t + c = — r 2 /3 
and t + c = — (r 2 + l/c)/3 where u r blows up and u vanishes. We can modify the solution 
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Figure 4. Graph of solution fl5 



by putting a cutoff on u at both parabolas so that 



u(t, r) 



t-p< -r 2 /3 or t - a > -r 2 /3 
r- 1/2 ^/~f(3(a - t) - r 2 )(3(t - J) + r 2 ), t - /3 > -r 2 /3 and t - a < -r 2 /3 



(see figure E]) where 



19) 



7 



3(a - /3) 



(90) 



and a = — c > (3 = — (c + l/3c) > with c > 0. Then the modified solution (189]) is C°° 
in r and t at all space-time points other than r = 0, a — t = r 2 /3 and /3 — t = r 2 /3. Near 

r = 0, for (3 < t < a, u = r -1 / 2 ^3a/7(« — t)(t — (3) + 0(r)j is singular (see figure [6]), while 
near the inner and outer cusps, u = (a/3(« — t) — r) 1 / 2 {^jQ^{a — (3) + 0(^3(a — t) — r) j 

for t < a and u = (r - y/3(P - t)) 1/2 ^67(0 - /3) + 0(r - ^/3(^ - t))) for t < /3 each 

have square-root behaviour in r (see figure [7]). From these properties, (155]) can be checked 
to satisfy the heat equation ( 176]) in a weak sense (i.e. u is only C° in t and r) on the spatial 
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domain r > 0. In particular, u has a finite, non-negative amount of heat 

'o. 



H 



t > a 



(3(a - t) - ry/ 2 (3(t - /3) + r^V^r 



T^A: 1/4 (3(a - /3)) 3/2 (3(« - t))- 1/4 2 Fi(l/4, 3/2; 3; {a - 0)/{a -t)), 0<t<a 
lb 



2 /2~ 



(3(a - t) - r 2 ) 1 / 2 (3(t - 0) + r^r^dr 



5 V vr 



r(3/4) V/ 4 (3(a - t)fl\F x {-l/2 : 3/2; 9/4; (a - t)/(a -/?)), t < (3 



(91) 



so thus u is in L 1 (IR + ) for all t > 0. The corresponding heat flux equation is given by 

0, t > a 

F — lim u/2 — F outCT + S 

r-s>0 

= -y A: 1/4 (3(a - /3)) 1 / 2 (3(a - t^^F^l/A, 1/2; 2; (a - f3)/(a - t)) 

- 2j P 1 (1/4, 3/2; 2; (a - /3)/(a - t))) , /3 < t < a 

-^inner -^outer ~t~ "S* 



!^r(3/4) 2 ^/4 (3(a - f))V*(5 aFi( _ 1/2) 3/2 . 5/4 . (Q _ t)/(a _ 0)) 

| 2^1(1/2, 3/2; 9/4; (a - f)/(a - /?))), 
a — p 7 



where 



F — lim -u/2 = lirn(rw r — u/2) = 0, 

-Pinner = 1™ TU r = +OO, 



outer 



s = < 



7-^3(0-*) 

l im r-u r 

r^A/3(Q-t) 

V 3 (°-*) 



-OO, 



a/ 3 ( q "*) 



(3(a - t) - r 2 )~ 3/2 (3(£ - 0) + r 2 )~ 3/2 r' 1/2 dr, (3 < t < 
(3(a - t) - r 2 )- 3 / 2 (3(t - 0) + r 2 )- 3/2 r~ 1/2 <*r, t < 0. 



a: 



t < (3 
(92) 

(93) 
(94) 

(95) 



(96) 



Here the quantities F innei . and F onteT are defined to be inward/outward heat fluxes arising 
from the inner and outer cusps, respectively. These fluxes act as cooling sources that cancel 
the singular contributions coming from the endpoints in the nonlinear heating source S in 
the flux equation (192]) . Their net effect gives a finite cooling rate, dH/dt < 0. Hence H 
decreases from its initial value H(0) ^ at t = to H(a) — at t — a. 
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Figure 6. Graph of solution (jS9J) with — 5, a = 20. 
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Figure 7. Graph of solution (J89J) with (3 = 5. 

As a result, non-similarity solution fl57|) physically describes an initial ring-shaped heat 
distribution in R 2 with radial temperature fronts at r = a/3(/3 — t) and r = a/3(o: — t) 
where < (3 < a. These fronts behave as heat flux sinks and move radially toward the 
origin at speeds dr/dt = — a/3/(/3 — t) and dr/dt = — y/3/(a — t). At time t = f3 the inner 
front reaches the origin and produces an infinite cooling point-source plus a compensating 
infinite outward heat flux that both persist until time t = a when the outer front reaches 
the origin. For all times < t < a, the heat absorbed at the temperature fronts exceeds the 
heat produced by the nonlinear source term u~ 3 at all points in R 2 inside the fronts, whereby 
the total amount of heat decreases to zero in a finite time t = a, with the heat distribution 
being given by u = at all points r > for times t > a. 

4. Concluding remarks 

As main results in this paper, analytically interesting exact solutions have been obtained 
for a mult i- dimensional semilinear heat equation ([1]) via a separation technique applied to 
the group foliation equations associated with the group of scaling symmetries ([2]) admitted 
by this equation. The solutions consist of explicit similarity solutions as well as other explicit 
solutions of a more general (non-similarity) form. 

In general our method provides a highly effective alternative to standard similarity re- 
duction for finding exact solutions to nonlinear PDEs with a group of scaling symmetries. 
Firstly, this method is an algorithmic refinement of the basic approach developed for the 
semilinear wave equation in [5], which leads to systematic reductions of the group foliation 
equations into overdetermined systems of ODEs that can be derived and solved by means 
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of computer algebra. In particular, these ODE systems are tractable to solve using the 
computer algebra package Crack [TT] . 

Secondly, the method is able to yield exact similarity solutions in an explicit form, whereas 
standard similarity reduction only gives an ODE that still has to be solved to find solutions 
explicitly and in general this step can be quite difficult. Indeed, the resulting similarity ODE 
02]) for the semilinear heat equation ([[]) cannot be solved by standard integration techniques 
such as symmetry reduction or integrating factors. 

Thirdly, because the group foliation equations contain all solutions of the given nonlinear 
PDE, our method can yield non-similarity solutions that are also not invariant under any 
other (non-scaling) point symmetries admitted by the nonlinear PDE. 

We can apply the same method more generally to nonlinear PDEs without scaling symme- 
tries by utilizing the group foliation equations associated with any admitted one-dimensional 
group of point symmetries of the given nonlinear PDE and by adapting the separation tech- 
nique to the specific form of the non- derivative terms that appear in the given group foliation 
equations. The algorithmic aspects of these steps will be the same as in the similarity case 
we have presented in this paper, since every one-dimensional group of point symmetries can 
be equivalently expressed as a group of scalings under an appropriate change of independent 
and dependent variables (i.e. by an invertible point transformation). 

For future work, we plan to present a full comparison between the present group foliation 
method and standard symmetry reduction as applied to many typical linear and nonlinear 
PDEs of interest, e.g. linear heat and wave equations; semilinear diffusion and telegraph 
equations; integrable semilinear evolution equations such as the Korteweg de Vries and 
Boussinesq equations. 
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